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SUMMARY 


Heat-transfer rates at the stagnation point of 1 .0- inch-diameter 
glass spherical models have been measured in shock-tube flows in air 
which corresponded to free-flight conditions between Mach numbers of 6.4 
and 15 -9 with stagnation temperatures up to 7^900° R. The heat-transfer 
rates were determined from measurements of the surface-temperatirre change 
with time of a thin-film-platinum resistance thermometer. The test results 
are presented and compared with the results from the theories of Lees 
(jet Propiilsion, April 1956) and Fay and Riddell (AVCO Research Report l) . 
The experimental results obtained give lower heat-transfer rates than 
both theories but are more compatible with results from the theory of Lees . 
These test results are also compared with thermocouple data and with the 
results obtained by the use of similar configurations. 


INTRODUCTION 


Many aerodynamic problems of hypersonic flight, such as the heat- 
transfer rate at the stagnation point, can be studied in a shock tube. 
This problem is of interest at present and is amenable to investigation 
because the shock tube can generate the high temperatirres necessary to 
simulate hypersonic conditions. For investigations of the heat-transfer 
rate at the stagnation point, only the stagnation enthalpy and pressure 
need be duplicated. The difference between the flight and test Mach 
numbers can be neglected if the external flow fields are similar. These 
conditions can be obtained in the straight portion of the tube for a 
limited flow region behind the primary shock wave . A series of tests 
were therefore made in a high-pressure shock-tube facility in the 
Langley Gas Dynamics Branch to obtain heat-transfer rates at the stag- 
nation point for a range of conditions available in a shock tube . These 
tests were made with a fast responding, thin-film-platinum resistance 
thermometer moiinted on the surface of a 1-inch-diameter glass model. The 
results of several tests are presented and show the heat-transfer rates 
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at the stagnation point obtained for conditions corresponding to a range 
in free-f light Mach number from 6.4 to 13.9 at high altitudes. The 
experimental results are compared with the theoretical results based on 
the work of Lees (ref. l) and Fay and Riddell (ref. 2). These results 
are discussed to show compatibility with the experimental results of 
Rose and Stark (ref. 3) and with the results of similar tests made with 
the use of a siurface thermocouple element of the same order of thickness 
as the thin-film thermometer. The thermocouple data were obtained in 
a previous investigation by Morton Cooper and Jim J. Jones of the 
Langley Aeronautical Laboratory. 


SYMBOLS 

c specific heat of glass, 0.20 Btu/lb-°F 

e current intensity, volts 

h enthalpy of air, Btu/lb 

i current, amp 

k diffusivity of glass, 6.67 x 10“^ sq. ft/sec 

M Mach number 

p press\ire, atm 

Q heat-transfer rate, Btu/sq ft/sec 

R resistance, ohms 

Reynolds number, based on model diameter 

t time, sec 

T temperature , 

T(t) surface-temperature rise as a function of time, 

p unit weight of glass, 139*2 Ib/cu ft 

Subscripts: 

0 stagnation region behind shock; also refers to sensitive element 

1 region ahead of primary shock 
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2 region behind primary shock 

5 region behind stationary model bow shock 

a assumed value 

e experiment 

s related to primary shock 

th related to theoiy 

w wall condition 


APPARATUS AND EXPERIMENTAL METHODS 

Heat-transfer measiirements were made in a 5f- - inch-diameter shock 

4 

tube with air as the test mediim. A shock in the tube was generated by 
bxirsting metal diaphragms with high-pressure heli\jm or hydrogen gas. 

When the resulting shock wave reached the model the temperature measure- 
ments began and continued for the duration of the approximately steady 
flow which followed. A more complete description of the shock tube, of 
the instrumentation, and of the flow details behind the primary shock 
wave is given in reference 4. The test models were located 68.15 feet 
from the diaphragm station and, in order to prevent shock reflections, 
the shock tube extended an additional 60 feet from the model. 

The wave velocity of the primary shock was measured at 10 locations 
in the tube by use of both ionization and pressirre sensitive probes. For 
conditions where weak shock waves failed to trip ionization probes, sev- 
eral ionization probes were replaced by pressvire transducers. Signals 
from all these instruments were detected by electronic interval counters 
or on an oscillograph to record time of shock passage. Shock strengths 
were varied by changing the driver gas, the driver-gas pressure, and 
the air pressure in the low-pressure section of the tube. Air pressures 

. » -P -2 

in the low-pressure section from approximately 0.45 X 10“ to 9*9 X 10 
atmosphere were used. 

Figure 1 is a photograph of typical models used. The models were 
formed into hollow 1.0- inch-diameter partial spheres at the end of a 
5/4-inch-outside-diameter and l/2-inch- inside-diameter pyrex glass tubing. 
In the vicinity of the stagnation point and extending rearward 45° from 
the stagnation point, the radius of the sphere had a tolerance of 
±0.01 inch. 
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The surface-temperature rise with time at the model stagnation 
point was determined hy measuring the resistance change of a sensitive 
element, or the thin-film thermometer. The sensitive element measured 
1/8 by 1/8 inch and consisted of a single thin coating of Hanovia Liquid 
Bright Platinum No. 05 which was applied to the model in the same manner 
as that described in reference 5* The resiolting sintered film had an 

approximate thickness of about 4 x 10“^ inch and an electrical resistance 
of about 16 ohms per square inch. This thermometer was resistant to 
abrasion by the flow and had a rapid response to temperatvure changes 
(ref. 6). Means were provided for conducting electrical signals from 
the sensitive element to instrumentation outside the shock tube through 
several leads that passed into the supporting tube of the model. (See 
fig. 1.) In order to minimize s\irface roughness behind the stagnation 
point, thin silver leads were applied to the surface to join the sensi- 
tive element with wire leads that entered the supporting tube . The sil- 
ver leads were made from a liquid silver paste and were applied with the 
same technique as that used for the platinum paste . The leads had a low 
resistance relative to the sensitive element. They were subsequently 
covered with one coating of clear insulating paint to help prevent 
detecting stray electrical signals. These insulated leads had an aver- 
age measured thickness of O.OOI5 inch. The wire leads in the supporting 
•tubing connected to the silver leads 90*^ behind the stagnation point, and 
four silver leads joined to the element. Two of these leads were current- 
carrying, while the other two were voltage-measuring leads - a technique 
used earlier in the investigation of reference 5 • Figiare 2 is a sche- 
matic drawing of the electrical circuit which connected to the sensitive 
element. In the circuit a 200-ohm resistor center tapped to ground was 
joined parallel to the thermometer to provide a balanced drain-off for 
any stray signal. 

The sensitive element was calibrated before each test by immersing 
the model in heated silicone oil and noting the resistance of the ele- 
ment for various temperatures of the oil. Dviring a test, the resistance 
change of the sensitive element on the model was measured by noting the 
voltage drop across the element for a current of approximately 10.0 milli- 
amperes through it. This current was found to produce \indetectable 
heating of the element when the model was placed in the testing environ- 
ment. The changing voltage at the sensitive element was amplified, fed 
into the vertical axis of a cathode-ray oscillograph, and recorded by a 
film-drimi camera. A camera film speed of 1.0 inch per millisecond was 
maintained by a synchronous motor. A typical voltage-time record pro- 
duced by the sensitive element is shown in figure 5 from which both the 
voltage rise Zse(t) and time elapsed t were obtained. 

As a result of bombardment of the model by fine particles in the 
flow, each model was used for only one test, although on occasions where 
the damage to the model was slight, as for the case of low-velocity flows. 


MCA TN 455A 


5 


the models were recalibrated for a second test. Fine particles were 
unavoidably present in the flow but it is believed that they reached 
the model only after the termination of steady flow. 

The thermocouple models of the previous tests of Morton Cooper and 
Jim J. Jones were solid glass spheres. The thermocouple was formed by 
evaporating onto the glass a compound band of overlapping nickel and 
silver which joined at the nose of the model. The electrical circuit 
joining the thermocouple with a voltage preamplifier and recorder con- 
sisted only of a series resistor which compensated for differences in 
resistance between the leads from the preamplifier to the ends of the 
thermocouple proper. The thickness of the thermocouple was roughly 

determined to be 4 x 10"^ inch. The method of calibrating the thermo- 
couple paralleled that of the thin-film thermometer. Unlike the thin- 
film thermometer, however, the thermocouple could be used for repeated 
tests without recalibration. The thermocouple readings were unaffected 
by abrasion and pitting during a test. All the tests with the theimio- 

couple model were made with an air pressure of O.566 x 10“^ atmosphere 
in the low-pressure chamber of the shock tube. 


DATA REDUCTION 

The surface -temperature rise at the stagnation point on the models 
was determined from the resistance change of the thin-film thermometer. 

If the relatively small resistance of the short wire leads to the instru- 
mentation is neglected the current through the sensitive element, before 
a test, is approximately 


1 - 1 

0 200 

where e represents the voltage observed at the sensitive element and 
i, the battery current measured by use of a standard resistor (fig. 2). 
The calibration of the sensitive element showed a linear resistance 
change with temperature and the test conditions did not exceed the upper 
calibrated temperature; therefore, the test-temperature change with time 
of the element could be obtained from the following equation: 


T(t) = Mil f^) 

^0 \^oj calibrated 

During a test the current iQ changed a negligible amount due to the 
resistance change of the sensitive element. 

In order to determine the heat-transfer rate at the stagnation 
point, a sensitive element of negligible thickness was used, that is. 
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approximately 4 x lO”^ inch (approximately 1.0 micron) . This thickness 
permitted the element to register nearly the true surface-temperature 
rise of the model (ref. 6). It also permitted, in computing the heat- 
transfer rate, approximation of the stagnation region as a semi-infinite 
solid when the running time was relatively short. With a large tempera- 
ture difference between the air and the model and with a region of approxi- 
mately steady flow behind the primary shock wave, the heat-transfer rate 
to the model for the first few instances of time could be considered con- 
stant . The heat-transfer rate could then be obtained from 

Q3 = ^ (1) 

2|ft 

As seen from equation (l) the behavior of the f\mction T(t) varies as 
the square root of time t . An actual voltage record is shown in figure 5 • 

Shock-tube flows generated by the primary shock are not ideally 
steady. Slight flow variations are sometimes present which affect the 
heat-transfer rate. Under these conditions, the heat-transfer rate can 
be described by 


Q(t) 


\/nkpc 

T(t) 1 

T(t) - t(t) 

Tt 


Jo (t-T)5/^ J 


( 2 ) 


which was derived from one-dimensional unsteady heat flow conditions. 
Although a constant heat-transfer rate was anticipated, the calculations 
were made by using eqmtion (2) as a verification. N\americal calcula- 
tion of this integral is complicated by the behavior of the integrand 
for values of the dumn^r variable t close to t. By terminating the 
integration at some value of t = tg^ < t, reasonable accuracy could be 
maintained when the rest of the integral was approximated by assuming 
that T was a linear function over the interval ta to t. This 
approximation leads to: 


Q(t) = 


i rtkpc 


T(t) 

\[t 


- 

2 Jo 


T(t) - T(t) 
(t - t) 


dT + 


dT(t) 


dt 


(t - ta) 


1/2 


(5) 


By comparing the altered integration method (eq. (5)) with the exact 
method for the case of constant Q,(t) (eq. (2)), a value of tg^ = 0.9t 

was found to maintain the accuracy of the nmerical integration but 
changed the value of Q,(t) by only k percent, a constant which could 
be accounted for in the result. Simpson's approximation rule was then 
used to compute the integral in equation (5), an<i the last term. 
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was the average slope between tg^ and t. Bulk properties of 

the particular glass were used for and c. The heat-transfer 

results obtained are presented in table l(a) with the calculated real- 
gas flow conditions and in figure 4 which shows the variation of heat- 
transfer rate with time. For comparison, table l(b) shows the flow con- 
ditions of air calculated for a constant ratio of specific heat of 1.4. 

As can be seen there is significant difference between the real and 
ideal stagnation-point conditions . 

Calculations of the flow conditions were based on the measirred pres- 
sure, temperature of the air ahead of the primary shock, and the shock 
velocity. Rankine-Hugoniot conservation relations were then used in an 
iteration process to determine all the flow conditions. Real-gas prop- 
erties for equilibrium conditions were used from reference 7 ‘to obtain 
table l(a) . The value of the dissociation energy of nitrogen used was 
incorrect but introduced only a small error in the calculated properties 
of air up to 7^900*^ R, since at the stagnation pressures the amount of 
nitrogen dissociation was small. 


The experimental values of the heat-transfer rate are compared in 
figure 5 with theoretical results from the theories of Lees (ref. l) and 
Fay and Riddell (ref. 2) . In the application of both theories a velocity 


gradient based on the Newtonian values 



was used, and real-gas 


properties were also considered in the theories. A Lewis number of 1.4 
was used in the theory of Fay and Riddell, and the theory of Lees was 


multiplied by 



RESULTS AND DISCUSSION 


After the primary shock wave passed the model, a region of compara- 
tively steady flow was established in which tests were made. The duration 
of this steady flow could be determined from the regularity of the voltage 
rise at the sensitive element. (See fig. 5*) This duration was longest 
for the low shock strengths, of the order of O.5 millisecond, and shortest 
for the high shock strengths, approximately 0.1 millisecond. Reference 4 
shows the running times available at other stations in the shock tube. 
Figure 4 shows the heat-transfer rates obtained at different time inter- 
vals . An apparent sharp rise in the heat-transfer rate, which is illus- 
trated as a typical case, was attributed to a particle hitting the model. 
The heat-transfer rates varied from 240 to 1,297 Btu/sq ft/sec for real- 
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gas stagnation temperatures that varied from 3^300° to 7^900° Since 
a longer test time tends to give more accuracy in the determined heat- 
transfer rate (see ref. 6) the heat rates at t = 0.22 millisecond were 
selected for comparison with theory for all hut the case of highest heat 
rate^ these rates are tabulated in table l(a) . Also included in this 
table are the results of calculations made by approximating the measured 
heat input as a constant rate . 

Figure 5 shows the comparison of the experimental heat-transfer 
rates obtained with those from application of the theories of Lees and 
Fay and Riddell by indicating the respective ratios. These ratios are 
plotted against free-flight Mach numbers which are the flight speeds 
required in a standard atmosphere (ref. 8) to duplicate the stagnation 
pressures and temperatures of the shock- tube tests . The test conditions 
corresponded to the Mach number range from M = 6.4 to 13-9 ^^nd to the 
altitude range between 40,000 and 95,000 feet. The maxim\mi heat-transfer 
rate of 1,297 Btu/sq ft/sec occurred with a stagnation condition that 
corresponded to a free-flight Mach nmber of 13 -9 • 

Figure 5 shows that the results obtained with the thin-film ther- 
mometer average from 0 to 30 percent less than the results calculated by 
use of the theory of Lees and from 20 to 40 percent less than the results 
calculated by use of the more accurate theory of Fay and Riddell. The 
thermocouple results of Cooper and Jones indicate better agreement with 
the theory of Fay and Riddell, the average of those resiilts being 10 per- 
cent less than the theoretical results. The res\ilts of reference 3 
agree favorably with the results from the theory of Fay and Riddell. 

There appears to be no apparent reason for the relatively lower values 
of the present test results . Precaution was taken to measure the heat- 
transfer rates at s\orface temperatures within the calibration range of 
the thin-film thermometer. As a note of interest, if the value of the 
Lewis number in the theory of Fay and Riddell were changed from 1.4 
to 1.0, this would lower the theory curve in figure 5 by not more than 

5 percent . 


CONCLUDING REMARKS 


Heat-transfer rates at the nose of 1.0-inch-diameter models were 
measured in a shock tube with flow stagnation temperatures up to 7,900° R 
by using a thin- film-platinum resistance thermometer. The maximum heat- 
transfer rate obtained was 1,297 Btu/sq ft/sec which corresponded to a 
free-flight Mach nirmber of 15 •9- 

The results of the present tests average approximately 20 to 40 per- 
cent lower than the results from the theory of Fay and Riddell (AVCO 
Research Report l) but are more compatible with the results from the 
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theory of Lees (jet Propulsion, April 1956), being from 0 to 30 percent 
lower. These test results also show lower heat-transfer rates than those 
obtained by use of a thermocouple and those obtained by Rose and Stark 
(AVCO Research Report 5) • 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics , 
Langley Field, Va., June l 8 , 1958* 
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TABLE I— TEST RESULTS 



(a) Real-gas conditions 


Test points 
from 
figure 4 

Pi» 

atm 

Tl, 

°R 

Ms 

P2» 

atm 

T2. 

°R 

«2 

%e,2 

P3. 

atm 

S’ 

°R 

M5 

P0> 

atm 

To. 

°R 

Btu/sq ft/sec 

Btu/sq ft/sec 
(a) 

0 

9.92 X 10"2 

558 

4.22 

2.08 

2,237 

1.71 

3.61 X 10^ 

6.70 

3,060 

0.628 

8.48 

3,236 

240 

226 

□ 

8.68 

510 

4.41 

1.99 

2,270 

1.73 

3.43 

6.56 

3,133 

.622 

8.25 

3,307 

307 

313 

0 

3.89 

539 

5.46 

1.39 

3,268 

1.92 

1.72 

5.62 

4,070 

.576 

6.90 

4,833 

543 

450 

A 

1.84 

540 

5.93 

.777 

3,716 

1.99 

.852 

3.40 

5,220 

.557 

4.00 

5,355 

487 

526 

V 

2.28 

532 

6.38 

1.12 

4,087 

2.05 

1.12 

5.22 

5,680 

.536 

6.20 

5,796 

491 

514 


.9^7 

536 

6.40 

.469 

4,120 

2.07 

.468 

2.26 

5,638 

.530 

2.63 

5,742 

477 

468 


1.99 

532 

6.69 

1.08 

4,362 

2.11 

1.02 

5.29 

5,980 

.520 

6.27 

6,109 

585 

581 

0 

2.04 

532 

7.24 

1.30 

4,842 

2.20 

1.13 

6.90 

6,512 

.499 

7.80 

6,622 

765 

684 

0 

1.01 

530 

7.35 

.633 

4,868 

2.25 

.579 

3.62 

6,445 

.487 

4.37 

6,570 

643 

625 

0 

1.99 

504 

7.84 

1.50 

5,135 

2.29 

1.26 

8.50 

6,880 

.489 

9.70 

6,998 

969 

862 

0 

.447 

523 

8.32 

.385 

5,364 

2.46 

.322 

2.48 

7,040 

.451 

2.78 

7,146 

692 

662 

a 

.618 

522 

9.14 

.646 

5,850 

2.58 

.508 

4.58 

7,843 

.440 

5.02 

7,934 

1,297 

1,356 


^Based on a constant rate of heat input. 


(b) Ideal-gas conditions, 7 a 1.4 


Test points 
from 
figure 4 

Pi. 

atm 

Tl. 

°R 

Ms 

P2. 

atm 

T2. 

OR 

M2 

%e,2 

P3, 

atm 

T3. 

°R 


Po. 

atm 

To. 

°R 

0 

9.92 X 10-2 

538 

4.22 

2.05 

2,368 

1.58 

3.18 X 10^ 

5.64 

3,260 

0.674 

7.65 

3,555 

□ 

8.68 

510 

4.41 

1.95 

2,404 

1.61 

3.02 

5.54 

3,345 

.667 

7.46 

3,642 

0 

3.89 

539 

5.46 

1.35 

3,627 

1.69 

1.36 

4.28 

5,274 

.642 

5.65 

5,708 

A 

1.84 

540 

5.93 

.752 

4,198 

1.72 

.651 

2.44 

6,156 

.635 

3.25 

5,685 

V 

2.28 

532 

6.38 

1.08 

4,711 

1.74 

.829 

3.64 

7,017 

.630 

4.75 

7,574 


.947 

536 

6.40 

.452 

4,781 

1.74 

.342 

1.53 

7,128 

.630 

2.00 

7,689 


1.99 

532 

6.09 

1.03 

5,130 

1.76 

.727 

3.54 

7,685 

.627 

4.63 

8,288 

Q 

2.04 

532 

7.24 

1.24 

5,917 

1.77 

.751 

4.35 

8,951 

.622 

5.65 

9,642 

0 

1.01 

530 

7.35 

.632 

6,003 

1.78 

.373 

2.22 

9,183 

.622 

2.89 

9,892 

0 

1.99 

504 

7.84 

1.43 

6,502 

1.79 

.782 

5.10 

9,913 

.619 

6.61 

10,675 

0 

.447 

523 

8.32 

.361 

7,536 

1.80 

.168 

1.31 

U,545 

.616 

1.69 

12,425 

A 

.618 

522 

9.14 

.601 

8,959 

1.82 

.231 

2.21 

13,824 

.615 

2.85 . 

14,870 
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Figure 1.- Shock- tube heat-transfer models. L-57-255-1 
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Figure 2.- Schematic diagram of electrical circuit used in measuring resistance change of 

stagnation-point thermometer. 
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Figure 5.- A typical record showing voltage rise with time at the sensi- 
tive element. Film speed, 1.0 inch per millisecond. 
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Duration of test, milliseconds 

Figure 4.- Observed heat-transfer rates for 0.5-inch nose radius. 
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I.- Heat-transfer rate at stagnation point compared with rates calculated from theory. 
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